Extended King Models for Star Clusters 

Kwok Sau Fa and I. T. Pedron 
Departamento de fi'sica, Universidade Estadual de Maringa, 
Av. Colombo, 5790, 87020-900 - Maringa - PR, Brazil 

Abstract 

Models applied to galaxies in equilibrium configuration are based on the 
solution of the collisionless Boltzmann equation and the Poisson equation for 
gravitational interaction which are related to each other by a smoothed-out 
mass density. This mass density is related to a distribution function /(x, v, t) 
and produces a gravitational potential (f>(x, t). The most popular models to 
describe spherically symmetric systems are King models. These models are 
based on a truncated isothermal sphere. They fit the surface brightness of 
globular clusters and some elliptical galaxies well. A well-known example 
of a galaxy that King models does not fit well is NGC 3379. In this work, 
we extend King models. Our models are based on the Tsallis distribution 
function. There is a parameter q in this function and we can recover the King 
distribution function in the limit q — > 1. We discuss the general behaviors of 
our models and, moreover, we use them to fit the surface brightness of NGC 
3379 and 47 TUC. 
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A lot of models for star clusters, based on the collisionless Boltzmann equation and 
Poisson equation, have been studied in the last century [l]j2|]. In this approach, the state of 
the system is described by a distribution function /(x, v,t) in phase space and considering 
the smooth gravitational potential. These models can be applied to star clusters in equilib- 
rium configuration or when the characteristic time of these structures does not exceed the 
relaxation time Jl|. 

The most commonly used models to describe spherically symmetric systems are King 
models. These models can fit the brightness distributions of globular clusters and some 
elliptical galaxies well |T|-§[. King models are based on the following distribution function 
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where e is the relative energy e = —E + $ = ^ — v 2 /2 and \I/ is the relative potential 
\[> = — $ + $o- This distribution function may be viewed as a truncated isothermal sphere. 
The characteristic of this model is that the relative potential \I/ and density of stars p become 
equal to zero at some finite radius r t (tidal radius). The value of r t depends on the central 



potential \l/(0)/o" 2 and King radius defined as r = y 9a 2 /AnGpo. In the limit ^(0)/a 2 —> oo, 
the isothermal sphere is recovered and r t it extends to infinity. 

In this work, we generalize King models. The main goal of this generalization is to try to 
describe some elliptical galaxies where King models cannot describe them well. Thus, we are 
also restricted to spherically symmetric models and the distribution function depends only 
on relative energy f(e). These kinds of models possess many solutions for the collisionless 
Boltzmann equation: Jeans stated that any steady-state solution can be expressed through 
integrals of motion on phase-space coordinates. Moreover, all these models are stable 

Our distribution function is based on Tsallis' statistics |J, with the entropy and dis- 
tribution function given by S = k(l — J2iPi)/(l ~ <?) an d = k (1 + (3(1 — g)e) , 
respectively. In this type of statistics, q represents a parameter which describes different 
statistics. It has been successfully applied to several systems such as Levi-type anomalous 
super diffusion |7], Euler turbulence f8j, and anomalous relaxation electron-phonon interac- 



tion [91. In astrophysics and cosmology we can mention the following set of works [ID 
but for our purposes, it is worthwhile to mention the work of Plastino and Plastino [10 
They have shown that Tsallis' entropy can determine the meaningful distribution function 
as described by polytropic models. Therefore, our distribution function, based on Tsallis' 
distribution function, is chosen as follows: 
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where q is a real parameter. This distribution function incorporates King's distribution 
function and recovers it in the limit q — ► 1. For simplicity, we also assume that the stars all 
have the same mass and the velocity distribution is isotropic everywhere. The dependence 
on r is determined by Poisson equation 

V 2 ^ = -AnGp , (3) 
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where p is the density given by 



p = An J f - -v 2 J v 2 dv . 

Substituting Eq. (f|) into (^)) we obtain 
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For convenience, we put the last equation in terms of \I//er 2 and dimensionless quantity 
r = r/ro, and it becomes 



A [r 2 — f *" 

dr \ dr Vct 2 , 



9r 2 [l + (1 - q) (J 



— M 



1-9 2 
IT 



d« - \ (£) 



/*(0) 



Jo 



*(0) 



(6) 



This equation should be integrated numerically. We use the Runge-Kutta method to calcu- 
late it. As usual, we choose the initial condition for <i\I/ /dr = at r = 0. Different models 
can be obtained by giving the central potential \E'(0)/cr 2 and the value of q. For q > 1, the 
validity of Eq. (^) is restricted by the condition (q — 1) < 1/ [^(0)/cr 2 ] in order to maintain 
the quantity inside of the square brackets positive. 

In order to compare the theory with observational data, the density calculated from Eqs. 
(|) and (H) is substituted into the projected density S(r) formula that is given by 
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where p(f) = p(f)/p - 

We should point out that King models are parameterized by the central potential ^(O) / a 2 
or the concentration c = log 10 (r t /r ). In these models, the tidal radius or c is uniquely 
determined by the value of the central potential. The greater the value of the central 
potential, the greater the tidal radius. In our models, the concentration is not uniquely 
determined by the central potential, but we can adjust the additional parameter q. In Fig.l 
we show the density profile for \l/(0)/<7 2 = 12 with different values of q. We can see that the 
tidal radius increases with the value of q. In Fig. 2 we show the surface brightness I(r) of 
our models, on the scale p = — 2. 5 log /(r) versus r 1 ' 4 , with different values of q. We also 
note that the curves become shallow with the decrease in the value of q, i.e., for q > 1, the 
curves become more crooked than those of q = 1, and for q < 1, the curves become more 
straight in some range and they become more and more straight with the decrease in the 
value of q. In the latter case, the larger the value of ^(0) / a 2 , the longer the straightend part 
of the curve becomes. This fact is important to the fitting of some elliptical galaxies, (for 
example NGC3379) because the surface brightness of King models cannot run straight over 
10 magnitudes with the increase in tidal radius and their curves become more crooked (see 
Fig. 2). In Fig.3, we fit the surface brightness of galaxy NGC3379. The value of r is chosen 
so that the theoretical curve falls on the observational data. We see that our theoretical 
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curve can fit the data of the outer region and runs almost straight over 10 magnitudes, in 
contrast to the King model which runs less than 10 magnitudes. In Fig. 4 we fit the surface 
brightness of globular cluster (47 TUC). We note that the King model (c=2.03) fits tolerably 
well, in contrast with our model which fits very well. 

In summary, we have extended King models to describe spherically symmetric systems. 
Our models are based on the distribution function (@) and can be viewed as a truncated 
Tsallis distribution function. In these models, we have considered that the stars all have 
the same mass and the velocity distribution is isotropic everywhere. We have analyzed the 
general properties of these models and contrasted them with King models (q = 1). As 
applications of our models , we have fitted the surface brightness of NGC3379 and 47 TUC. 
From Fig. 3, we see that the model with g ^ 1 fits the surface brightness of the galaxy 
tolerably well, in contrast to the case of q = 1 which cannot fit the outer region. Moreover, 
the improvement of our model in relation to the King model is clear from Fig. 4. These 
results suggest that the velocity distribution of stellar systems described here deviates from 
Maxwellian shape. Of course, our models represent the simplest of models, but can be 
extended to include the different masses and anisotropy of actual clusters. 
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FIGURE CAPTIONS 



Fig.l - Density profile plots with \l/(0)/cr 2 = 12 and different values of q. 

Fig. 2 - Surface brightness plots versus r 1//4 with \I/(0) / cr 2 = 8, 12 and different values of q. 
The dotted curves correspond to \I/(0)/(X 2 = 8: from bottom to top the value of q is given by 
q = 0.92, 0.96, 1, 1.02, 1.04. While, the full curves correspond to^(0)/(T 2 = 12: from bottom 
to top the value of q is given by q = 0.6, 0.8, 1, 1.02, 1.04. 

Fig. 3 - Surface brightness fit of galaxy NGC3379 measured by different observers along 
the galaxy's east-west axis (from data published by Vaucouleurs and Capaccioli f§]). The 
dashed curve corresponds to \l/(0)/cr 2 = 9.2 and q — 1, whereas the full curve corresponds 
to ^(0)/a 2 = 12.2 and q = 0.939. 

Fig. 4 - Surface brightness fit of globular cluster 47 TUC by different observers (from 
data published by Illingworth & Illingworth Q and King et. a!., 1968 Q). The dotted 
curve corresponds to \l/(0)/er 2 = 8.868 and q = 1, whereas the full curve corresponds to 
^(0)/a 2 = 10 and q = 0.966. 
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